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Abstract. A Hamiltonian defined as a polynomial in creation and annihi- 
lation operators is considered. After a minimization of its expectation value 
over pure Gaussian states, the Hamiltonian is Wick-ordered in creation and 
annihillation operators adapted to the minimizing state. It is shown that 
this procedure eliminates from the Hamiltonian terms of degree 1 and 2 
that do not preserve the particle number, and leaves only the terms that 
can be interpreted as quasiparticles excitations. We propose to call this 
fact Beliaev's Theorem, since to our knowledge it was mentioned for the 
first time in a paper by Beliaev from 1959. 



1. Introduction 

Various phenomena in many-body quantum physics are explained with help 
of quasiparticles. Unfortunately, we are not aware of a rigorous definition of 
this concept, except for some very special cases. 

A typical situation when one speaks about quasiparticles seems to be the 
following: Suppose that the Hamiltonian of a system can be written as if = 
Hq + V, where Hq is in some sense dominant and is a perturbation that in 
first approximation can be neglected. Suppose also that 



where i? is a number, operators b*/bi satisfy the standard canonical commuta- 
tion/ ant icommut at ion relations (CCR/CAR) and the Hilbert space contains 
a state annihilated by bi (the Fock vacuum for b.i). We then say that the 
operators b*/bi create /annihilate a quasiparticle. 
Of course, the above definition is very vague. 

In our paper we describe a simple theorem that for many Hamiltonians 
gives a natural decomposition H = Hq + V with Hq of the form (11. ip . and 
thus suggests a possible definition of a quasiparticle. Our starting point is a 
fairly general Hamiltonian H defined on a bosonic or fermionic Fock space. 
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For simplicity we assume that the 1-particle space is finite dimensional. With 
some technical assumptions, the whole picture should be easy to generalize to 
the infinite dimensional case. We assume that the Hamiltonian is a polynomial 
in creation and annihilation operators a*/ai, i = 1, . . . ,n. (This is a typical 
assumption in Many Body Quantum Physics and Quantum Field Theory). 

An important role in Many Body Quantum Physics is played by the so-called 
Gaussian states, called also quasi-free states. Gaussian states can be pure or 
mixed. The former are typical for the zero temperature, whereas the latter for 
positive temperatures. In our paper we do not consider mixed Gaussian states. 

Pure Gaussian states are obtained by applying Bogoliubov transformations 
to the Fock vacuum state (given by the vector Q annihilated by Oj's). Pure 
Gaussian states are especially convenient for computations. 

We minimize the expectation value of the Hamiltonian H with respect to 
pure Gaussian states, obtaining a state given by a vector Cl. By applying an 
appropriate Bogoliubov transformation, we can replace the old creation and 
annihilation operators a*, ai by new ones 6*, 6j, which are adapted to the "new 
vacuum" fi, i.e., that satisfy hiVl = 0. We can rewrite the Hamiltonian H in 
the new operators and Wick order them, that is, put b* on the left and bi on 
the right. The theorem that we prove says that 

H = B + J2 Dijb*bj + V, 

ij 

where V has only terms of the order greater than 2. In particular, H does not 
contain terms of the type b*, bi, b*b*, or bibj. It is thus natural to set Hq := 
B + 'Yliij Dijb*bj. Dij is a hermitian matrix. Clearly, it can be diagonalized, so 
that Hq acquires the form of (11. ip . 

We present several versions of this theorem. First we assume that the Hamil- 
tonian is even. In this case it is natural to restrict the minimization to even 
pure Gaussian states. In the fermionic case, we can also minimize over odd 
pure Gaussian states. In the bosonic case, we consider also Hamiltonians with- 
out the evenness assumption, and then we minimize with respect to all pure 
Gaussian states. 

The procedure of minimizing over Gaussian states is widely applied in prac- 
tical computations and is known under many names. In the fermionic case in 
the contex of nuclear physics it often goes under the name of the Hartree-Fock- 
Bogoliuhov method [11]. It is closely related to the Bardeen- Cooper- Schrieffer 
approximation used in superconductivity [1] and the Fermi liquid theory devel- 
oped by Landau [lOJ- In the bosonic case it is closely related to the Bogoliubov 
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approximation used in the theory of superfluidity [1] , see also fl2[ [5] . In both 
bosonic and fermionic cases it is often called the mean- field approach [S]. 

The fact that we describe in our paper is probably very well known, at 
least on the intuitive level, to many physicists, especially in condensed matter 
theory. One can probably say that it summarizes in abstract terms one of the 
most widely used methods of contemporary quantum physics. The earliest 
reference that we know to a statement similar to our main result is formulated 
in a paper of Beliaev . Beliaev studied fairly general femionic Hamiltonians 
by what we would nowadays call the Hartree-Fock-Bogoliubov approximation. 
In a footnote on page 10 he writes: 

The condition H20 = may be easily shown to be exactly equivalent to the 
requirement of a minimum "vacuum" energy U . Therefore, the ground state 
of the system in terms of new particles is a "vacuum" state. The excited states 
are characterized by definite numbers of new particles, elementary excitations. 
Therefore, we propose to call the main result of our paper Beliaev 's Theorem. 

The proof of Beliaev's Theorem is not difficult, especially when it is is for- 
mulated in an abstract way, as we do. Nevertheless, in concrete situations, 
when similar computations are performed, consequences of this result may of- 
ten appear somewhat miraculous. The authors of this work witnessed it several 
times: the authors themselves, or their colleagues, after tedious computations 
and numerous mistakes watched the unwanted terms disappear |5l [6j . As we 
show, these terms have to disappear by a general argument. 
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2. Preliminaries 

2.1. 2nd quantization. We will consider in parallel the bosonic and fermionic 
case. 

Let us describe our notation concerning the 2nd quantization. We will al- 
ways assume that the 1-particle space is C". (It is easy to extend our analysis 
to the infinite dimensional case). The bosonic Fock space will be denoted 
rs(C") and the fermionic Fock space ra(C"). We use the notation Ts/ai'C"') for 
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either the bosonic or fermionic Fock space. Q G rs/a(C") stands for the Fock 
vacuum. If r is an operator on C", then r(r) stands for its 2nd quantization, 
that is 



55" 



r(r) := © r 

a*, ai denote the standard creation and annihilation operators on Fs/alC"), 
satisfying the usual canonical commutation/anticommutation relations. 

2.2. Wick quantization. Consider an arbitrary polynomial on C", that is a 
function of the form 



(2.2) h{-z,z):=Y,hc 

where 2; = (zi, . . . , z„) G C", z denotes the complex conjugate of z and a = 
(«!,..., a;„) G (N U {0})" represent multiindices. In the bosonic/fermionic 
case we always assume that the coefficients ha^p are symmetric/antisymmetric 
separately in the indices of z and z. 

We write |a| = ai + ■ ■ ■ + a;„. We say that h is even if the sum in (12. 2 p is 
restricted to even |a| + . 

The Wick quantization of (12. 2 p is the operator on rs/a(C") defined as 

(2.3) hia*,a) := J^haA^^T^^ ■ 

In the fermionic case, (12. 3 p defines a bounded operator on ra(C"). In the 
bosonic case, (12. 3 p can be viewed as an operator on f] DomA^" C rs(C"), 

n>0 

where 



a- ai 



i=l 



is the number operator. 



2.3. Bogoliubov transformations. We will now present some basic well 
known facts about Bogoliubov transformations. For proofs and additional 
information we refer to |3] (see also [7], [9]). We will often use the summation 
convention of summing with respect to repeated indices. 
Operators of the form 

(2.4) Q = 6ija*a* + h^ialai + Oijajai, 

where /i is a self-adjoint matrix, will be called quadratic Hamiltonians. In the 
bosonic/fermionic case we can always assume that 9 is symmetric/antisymmetric. 
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The group generated by operators of the form e^'^, where Q is a quadratic 
Hamiltonian, is called the metaplectic (Mp) group in the bosonic case and the 
Spin group in the fermionic case. 

In the bosonic case, the group generated by Mp together with e^^^'"**"'"^''^'^ 
Hi & C, i = 1, . . . ,n, is called the ajfine mataplectic (AMp) group. 

In the fermionic case, the goup generated by operators yia* + y^ai with 
^ = 1 (which are unitary) is called the Pin group. Note that Spin is a 
subgroup of Pin of index 2. 

In the bosonic case, consider U G AMp. It is well known that 

(2.5) U UiU* = pijaj + Qija* + ^i, Ua*U* = p^^a* + q^jaj + 

for some matrices p and q and a vector ^. 

In the fermionic case, consider U G Pin. Then 

(2.6) UaJJ* = pijttj + qijtt* , Ua*U* = p^^a* + q^^aj 

for some matrices p and q. 

The maps ( 12. Sp and ( 12. 6p are often called Bogoliubov transformations. Bo- 
goliubov transformations can be interpreted as automorphism of the corre- 
sponding classical phase space. Let us describe briefly this interpretation. 

Consider the space C" © C". It has a distinguished 2?T,-dimensional real 
subspace consisting of vectors {z,z) = ((-2i)j=i,...,n, (^i)j=i,...,n), which we will 
call the real part of C" © C", and which can be interpreted as the classical 
phase space. The real part of C" © C" is equipped with a symplectic form 

(2.7) {z,z)uj{z',z') ■.= lm{z\z'), 
and a scalar product 

(2.8) {z,z) ■ {z',z') ■.= Re{z\z'). 

Consider the bosonic case. Note that the transformation (12. 5p . viewed as a 
map on the real part of C" © C" given by the matrix 



p q 
q p 



and the vector 

, preserves the symplectic form (12. 7p - in other words, it belongs to ASp, 



the affine symplectic group. More precisely, it is easily checked that in this 
way we obtain a 2-fold covering homomorphism of AMp onto ASp. 

In the fermionic case there is an analogous situation. The transformation 

p q 



(12.61) . viewed as a map on the real part of C"©C" given by the matrix 



q p 
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preserves the scalar product (12 .Sp - in other words, it belongs to O, the or- 
thogonal group. More precisely, it is easily checked that in this way we obtain 
a 2-fold covering homomorphism of Pin onto O. 

2.4. Pure Gaussian states. We will use the term pure state to denote a 
normalized vector modulo a phase factor. In particular, we will distinguish 
between a pure state and its vector representative. 

On Fock spaces we have a distinguished pure state called the (Fock) vacuum 
state, corresponding to Q. States given by vectors of the form UQ, where 
U G Mp OT U E Spin, will be called even pure Gaussian states. The family of 
even pure Gaussian states will be denoted by C5s/a,o- 

In the bosonic case, states given by vectors of the form UQ where U G AMp 
will be called Gaussian pure states. The family of bosonic pure Gaussian states 
will be denoted by 6s- 

In the fermionic case, states given by vectors of the form UQ, where U G Pin 
will be called fermionic pure Gaussian states. The family of fermionic pure 
Gaussian states is denoted 6a- 

Fermionic pure Gaussian states that are not even will be called odd fermionic 
pure Gaussian states. The family of odd fermionic pure Gaussian states is 
denoted 0a,i- 

One can ask whether pure Gaussian states have natural vector representa- 
tives (that is, whether one can naturally fix the phase factor of their vector 
representatives). In the bosonic case this is indeed always possible. If c = [cij] 
is a symmetric matrix satisfying ||c|| < 1, then the vector 

(2.9) det(l - c*c)^/^e^'=»^<"Jn 

defines a state in C5s,o (see [13j). If = [Oij] is a symmetric matrix satisfying 



tanh ' 



^-e, then (^M> equals 



c = 1 

(2.10) e'^m 
with 

(2.11) Xg := 6ija*a* + Oijajai. 

Each state in 0s,o is represented uniquely as (12. 9 p (or equivalently as (I2.10p ). 
In particular, (I2.10p provides a smooth parametrization of (Ssfl by symmetric 
matrices. 

The manifold of fermionic even pure Gaussian states is more complicated. 
We will say that a fermionic even pure Gaussian state given by \E' is nonde- 
generate if (^^I^P) 7^ (if it has a nonzero overlap with the vacuum). Every 
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nondegenerate fermionic even pure Gaussian state can be represented by a 
vector 

(2.12) det(l + c*c)-^/'e^"'^<"^*fi, 

where c = [cjj] is an antisymmetric matrix, li 6 = [6ij] is an antisymmetric 
matrix satisfying c = i^-S^0^9, \\9\\ < n/2, then (1212]) equals 

(2.13) e'^'n 
with 

(2.14) Xg := OijOla* + ^a^ai. 

Not all even fermionic pure Gaussian states are nondegenerate. Slater de- 
terminants with an even nonzero number of particles are examples of even 
Gaussian pure states that are not nondegenerate. Note that vectors fl2.12p are 
natural representatives of their states. It is easy to see that only nondegenerate 
fermionic pure Gaussian states possess natural vector representatives. 

Nondegenerate pure Gaussian states form an open dense subset of 0a,o con- 
taining the Fock state (corresponding to c = 6* = 0). In particular, fl2.12p 
provides a smooth parametrization of a neighborhood of the Fock state in 0a,o 
by antisymmetric matrices. 

The fact that each even bosonic/nondegenerate fermionic pure Gaussian 
state can be represented by a vector of the form (12. 9p / (I2.12p goes under the 
name of the Thouless Theorem. (See [H]; this name is used eg. in the mono- 
graph by Ring and Schuck [11]). The closely related fact saying that these 
vectors can be represented in the form (12.101) / (12.131) is sometimes called the 
Ring-Schuck Theorem. 

By definition, the group AMp/Pin acts transitively on ©s/a- In other words, 
for any Cl G ©s/a we can find U G AMp/Pin such that Cl = UQ. Such a. U is 
not defined uniquely - it can be replaced by Ur{r), where r is unitary on C". 

Clearly, if we set 

(2.15) bi := UaiU*, b* := Ua*U\ 

then biVL = 0, 2 = 1, . . . , n, and they satisfy the same CCR/CAR as a,, i = 
1, . . . ,n. If /i is a polynomial of the form (12. 2p . then we can Wick quantize it 
using the transformed operators: 

hib*,b) = J2haAbTb^- 

a, 13 

Obviously, Uh{a*,a)U* = h{b*,b). 
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3. Main result 

As explained in the introduction, we think that the following result should 
be called Beliaev's Theorem. 

Theorem 3.1. Let h be a polynomial on C" and H := h{a*,a) its Wick 
quantization. We consider the following functions: 

(1) (bosonic case, even pure Gaussian states) 0s,o 9 $ 

(2) (bosonic case, arbitrary pure Gaussian states) (Sg 3 ^ ^ 

(3) (fermionic case, even pure Gaussian states) C5a,o 9 

(4) (fermionic case, odd pure Gaussian states) (Sa,i 9 $ i-^ 

In (1), (3) and (4) we assume in addition that the polynomial h is even. 

For a vector Cl representing a pure Gaussian state, let U e AMp/Pin satisfy 
Vl = Ufl. Set bi = UaiU* and suppose that h is the polynomial satisfying 
H = h{b*, b). Then the following statements are equivalent: 

(A) represents a stationary point of the function defined in (l)-(4). 

(B) 

h{b* , b) — B + Dijb*bj + terms of higher order in 6's. 

Proof. Let us prove the case (2), which is a little more complicated than the 
remaining cases. Let us fix ?7 G AMp so that Cl = UQ. Clearly, we can write 

H = h{b\ b)^B + Kibi + Kib* + Oijb]b* + Oijhb, + 

(3.16) +terms of higher order in 6's. 

We know that in a neighborhood of Vt arbitrary pure Gaussian states are 
parametrized by a symmetric matrix 9 and a vector y: 

^ C/e^^(2')e^^«Q, 

where Xg :— 6ija*a* + OijUjai and (f){y) — y^a* + y^a^. We get 

{Ue"t'^y^fi'^on\HU(i"t'^y^e'^''Q) = {e"*'^y^e'^'>n\U*h{b\ b)Ue"^^y^e'^'>n) 

(3.17) = (fi|e-'^»e-*'^(J')/i(a*,a)e'*('^)e*^«fi). 
Now 

e-i^«e-'^(^)^(a*,a)e^<^(^)e'^'' = S - i(%Oy - %0i,) - i(^,K, - t/,Z,) 

+terms containing Oj or a* + 0(||^||^, \\y\\^). 
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Therefore, (13.171) equals 

(3.18) B - i(%0,, - %0,,) - i{y,K, - yj{,) + 0{\\ef, \\y\\^). 

Since vectors y and matrices 6 are independent variables, (I3.18P is stationary 
at fl if and only if [Oij] is a zero matrix and [Ki] is a zero vector. This ends 
the proof of part (2). 

To prove (3) and (4) we note that for U G Pin, the neighborhood of Q = UQ 
in the set of fermionic pure Gaussian states is parametrized by antisymmetric 
matrices 6: 

where again := 6ija*a* + Oijajtti. Therefore, it suffices to repeat the above 
proof with yi = Ki = 0, i = 1, . . . ,n. 
The proof of (1) is similar. 

Proposition 3.2. In addition to the assumptions of Theorem \3.1\ (2), suppose 
that Vt corresponds to a minimum. Then the matrix [Dij\ is positive. 

Proof. Using that O and K are zero, we obtain 
e-^<^(2')/i(a*,a)e''^(2/) = B + %Dijyj 

+terms containing aj or a* + 0(||2/||^). 

Therefore, (13.171) equals 

(3.19) B + y^D,,y, + 0{\\yf). 

Hence the matrix [Dij] is positive. 

Note that in cases (1), (3) and (4) the matrix [Dij] does not have to be 
positive. 
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